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Abstract 

Charged lepton transverse momenta in the Drell- Yan processes play an impor- 
tant role at the LHC in precision measurements of the Standard Model parameters, 
such as the VF-boson mass and width, their charge asymmetries and sin^ 9w- There- 
fore, their distributions should be described as accurate as possible by the Monte 
Carlo event generators. In this paper we discuss the problem of matching the hard- 
process kinematics of the Monte Carlo generator WIN HAG with the parton-shower 
kinematics of the PYTHIA 6.4 generator while interfacing these two programs. We 
show that improper assignment of the quark and antiquark effective momenta in 
the LO matrix element computations may affect considerably the predicted lep- 
ton transverse momenta and even completely reverse their charge asymmetries at 
the LHC. We propose two matching schemes in which the NLO QCD distributions 
of the leptonic kinematical variables can be well reproduced by the LO WIN HAG 
generator. 
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1 Introduction 



Distributions of transverse momenta of charged leptons (pip) produced in Drell-Yan (DY) 
processes are important observables in hadron collider experiments. Their sensitivity to 
the values of the Standard Model (SM) parameters and to the polarisation of W and Z 
bosons produced in the DY process can be used in precision tests of the Standard Model 
(SM) and in searches for new physics. 

The measurement of the SM parameters: VT-boson mass (Mw), its charge asymmetry 
{Myy+ — M\Y-) and width iVw) at the LHC was investigated in Refs. (l}]6]. It was shown 
that in order to reach a precision at the level of ~ 10 MeV, the p^rp distributions must be 
controlled experimentally to a comparable accuracy. Since the Monte Carlo (MC) event 
generators are indispensable tools to derive the values of the SM parameters from the 
measured distributions, they must predict these distributions to even higher precision. 

MC event generators for Drell-Yan processes developed so far can be divided into two 
categories. The first one includes generators which are based on the precise calculations 
of matrix elements (ME) for hard processes: those including QCD effects to the NLO 
or even NNLO level, and those including QED and electroweak (EW) corrections. In 
these generators, called hereafter the ME generators, differential cross sections for the 
hard process are convoluted with the universal parton distribution functions (PDFs) of 
the hadron beams. These PDFs depend, apart from the factorisation scale, only on 
the longitudinal momenta of partons, x. In this type of the MC generators partons 
entering the hard process are assumed to be parallel to hadron beam direction. The 
second category of MC generators includes the so-called parton-shower (PS) generators, 
such as PYTHIA [T], HERWIG [8], etc. They generate initial-state multi-parton radiation in 
form of the LO-type QCD/QED parton cascade and then perform hadronisation as well as 
some particle decays. In the PS generators, partons entering the hard process are no longer 
coUinear with the hadron-beam directions, but acquire non-zero transverse momenta. The 
hard process itself is described by these generators usually at the LO level. Thus, as long 
as the ME precision is the key factor determining the overall accuracy of the measurement 
of a selected observable, they are inferior with respect to the ME generators. However, 
in the remaining cases they often provide better description of the hadronic energy flow 
associated with the DY process, in particular for not to high transverse momenta of the 
W and Z bosons. 

A very important, and at the same time difficult issue is how to combine these two 
types of generators, avoiding, on the one hand, a double counting of QCD corrections 
and, on the other, possible gaps in phase space present in some PS algorithms. As the 
state-of-the-art practical solutions to this problem for the QCD effects are regarded the 
MC@NLO |9| and POWHEG (lO) generators. They match the NLO QCD ME calculations 
with the parton shower generators, albeit in different ways. MCONLO uses HERWIG 



or HERWIG+-I- 11 for parton shower generation, while POWHEG is more universal, in 



principle it can use an arbitrary parton-shower generator. 

In the case of ME generators that include QED/EW radiative corrections the situation 
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is simpler as combining them with PS generators usually does not lead to double counting 
One of such generators is the MC program WINHAC 12 14 . It includes higher-order 
QED effects for the final-state radiation (FSR) and initial-final state interferences in the 
Yennie-Frautschi-Suura (YFS) exclusive exponentiation scheme 15 together with the 
0{a) EW corrections for the full charged- current DY process 16 . For the QCD (and 
QED ISR) effects it is interfaced with PYTHIA 6.4 ff]. However, we do not use the Les 
Houches Accord (LHA) scheme y/7j but our own interface, in which the PYTHIA routines 
for parton-shower generation and hadronisation are called directly from the WINHAC 
program. The principal reason is that the LHA is not general enough to transmit the full 
information contained in the spin density matrix of the W and Z bosons between these 
two generators. Moreover, from the purely technical perspective, we avoid writing/reading 
events into/from disk files, which makes event generation easier and much more efficient. 
For example, in our studies presented in Refs. [l]-[6] and requiring generation of C(10^^) 
events, its efficiency was one of the principal optimisation targets. 

In any interface which extends the LHA scheme to processes involving the spin-1 EW 
bosons as intermediate particles a particular care must be taken of the proper matching of 
the ME-type kinematics with the PS-type kinematics. In the following we show that this is 
particularly important to describe the charged lepton transverse momenta distributions in 
the Drell-Yan processes at the LHC In particular, we find that using the original PYTHIA 
effective momenta of incoming quarks in the above matching results in strongly biased 
pip distributions, which is particularly visible in their charge asymmetries. We propose 
solutions to this problem that seem to reproduce well the NLO QCD predictions for these 
asymmetries, as obtained e.g. from MCONLO. The goal of this exercise is to try to achieve 
the NLO QCD precision for the description of the leptonic observables in the DY process 
with the suitably matched LO QCD generator which incorporates the state-of-the-art EW 
corrections. 

The paper is organised as follows. In the next section we describe in detail the matching 
of the ME kinematics with PS kinematics as it is realised in the WINHAC MC event 
generator. In Section 3 we present numerical results illustrating the above issues, discuss 
their meaning for the LHC physics and propose our solutions that match the NLO QCD 
predictions. Finally, in Section 4 we summarise the paper. 



2 Matching of ME and PS kinematics 

The WINHAC |12j MC event generator is dedicated to precision modelling of single W- 
boson production with leptonic decays, i.e. the charged-current Drell-Yan processes, in 
proton-proton, proton-antiproton and ion-ion collisions, with the main emphasis on the 
QED effects and electroweak corrections. It uses fully massive spin amplitudes to evaluate 
the hard process matrix elements. They can be computed in an arbitrary reference frame, 

^The only problem here may be the QED ISR, but since its numerical effects are rather small, it can 
be dealt sufficiently well by PS generators. 
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in particular they can be used to calculate polarised H^-boson cross section^ In terms of 
the perturbative QCD the current version includes the LO hard-process matrix element. 
The QCD effects enter only through scaling-violating PDFs, taken from the LHAPDF 
library [Ts]. Therefore, in WINHAC incoming quarks producing a VT-boson are coUinear 
with hadron (ion) beams; their longitudinal momenta are given by the two x-variables 
which are generated according two PDFs and subsequently convoluted with the hard- 
process differential cross section. At this stage WINHAC is a ME-type MC generator. 
Its full event kinematics (i.e. all four-momenta of initial, intermediate and final state 
particles) is constructed for incoming partons coUinear with the colliding beams. Let us 
call it the ME kinematics. At this level WINHAC has been cross-checked numerically to a 



high precision with independent calculations 13 ,16 ,19 . 

Events of the ME kinematics do not look very realistic from the experimental point 
of view for the following two reasons. Firstly, QCD radiation affect not only longitudinal 
momenta of partons but also their transverse momenta. Therefore, using purely coUinear 
PDFs for the description of the QCD effects is not sufficient. Secondly, partons are not 
observed experimentally. What can be observed are the products of their hadronisation 
and decays. Therefore, in a realistic MC generator, to be used in an experimental data 
analysis, the above effects must be taken into account. In WINHAC this is done through 
the interface to the PYTHIA 6.4 generator which performs the initial-state LO-type QCD 
(and QED) parton shower, appropriate proton-remnant treatment, and necessary hadro- 
nisation/decays. In PYTHIA, partons entering the hard process are not coUinear with the 
hadron beams. In the case of the charged-current DY processes, PYTHIA provides, in its 
event record, the momenta of the two effective on-mass-shell quarks producing VT-bosons. 
A vectorial sum of these momenta gives a momentum of an appropriate VT-boson. Such 
a ly-boson, in contrast to the case of a ME-type MC generator, carries the transverse 
momentum, being a vectorial sum of the quarks transverse momenta. This has to be 
taken into account in constructing the hard-process event kinematics. Let us call the 
kinematics in which the incoming partons are not parallel to the beams (as a result of the 
aforementioned effects) the PS kinematics, and a corresponding MC generator - the PS 
generator. 

While interfacing the ME-type generator with the PS generator one has to take care 
of the appropriate matching of the ME kinematics with the PS kinematics. This is par- 
ticularly important for processes in which particles with non-zero spin, e.g. IV/Z-bosons, 
are produced as intermediate states. In the following, we describe in detail how such a 
matching is performed in the WINHAC interface to PYTHIA 6.4. Then we discuss possible 
pitfalls of the kinematical matching of these two types of MC generators. 

In the WINHAC interface to PYTHIA the final hard-process event kinematics is con- 
structed through the foUowing steps: 

1. The ME kinematics of a given MC event is generated in WINHAC- the four-momenta 
of: the incoming quarks, the intermediate VT-boson, the final-state leptons, and the 



^WINHAC provides several options for computing polarised W^-boson processes in several reference 
frames. 
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radiative photons are constructed in the frame in which quarks are coUinear with 
the hadron beams, where the +z axis is the direction of one of the beams at the 
colhsion point. 

2. All the above four-momenta are then Lorentz-boosted along the W-hoson direction 
to the ly-boson rest frame, which is also the centre-of-mass frame of the incoming 
quarks. In this frame quarks are still aligned along the z axis. 

3. The PS kinematics is generated by PYTHIA, in which the effective on-mass-shell 
quarks producing the W-bo^on are non-collinear with the beams as a results of the 
initial-state QCD/QED parton showei|^ Their four- momenta are given in the LAB 
frame with the +z axis along one of the hadron beams. 

4. The above PYTHIA-quarks four- momenta are Lorentz-boosted along the sum of 
their momenta to their centre-of-mass frame, which is also the VT-boson rest fram^ 
Contrary to the WINHAC quarks in point 2, their momenta, although back-to-back, 
are not aligned with the z axis in this frame. Instead, their direction is rotated with 
respect to the z axis by the polar angle Oq, and with respect to the x axis by the 
azimuthal angle 0^. 

5. From the PYTHIA-quarks momenta, specified in the above frame, we calculate the 
angles {9q, (pg) and then we perform rotations of all the WINHAC momenta specified 
in the H^-rest frame (point 2) using these rotation angles. After such rotations the 
WINHAC quarks are aligned with the PYTHIA quarks. 

6. Finally, the whole WINHAC event is Lorentz-boosted from the above frame to the 
LAB frame along the sum of the original PYTHIA-quarks momenta (the boost is 
opposite to the one in point 4). 

In our opinion all the above steps are needed for a proper matching of the ME kinematics 
with the PS kinematics in any interface between the ME-type MC event generator and the 
PS generator. This is particularly important for the production of the VF-bosons which 
are spin-1 particles with V — A couplings to fermions. In such a case respecting all the 
spin correlations in the above matching is obligatory. 

The above kinematical matching relies on the correct PYTHIA generation of the in- 
coming "effective" quarks momenta. Their spacial orientation is crucial for the spin 
correlations which, in turn, influence the angular distributions of the W-decay leptons 
and, as a consequence, their distribution. 

At this point one may ask if such effective on-shell quarks make sense at all. It 
has been known for some time that a cross section corresponding to the real NLO QCD 

•^Actually, for technical reasons PYTHIA performs the so-cahed backward QCD evolutfon. This aspect 
is not important for our discussion. Here we are concerned mainly with the PS kinematics in which 
partons entering the hard process are not parallel to the hadron beams - this can be a result of any type 
of a parton-shower algorithm in which transverse degrees of freedom are not neglected (integrated out). 

^Instead of a single parallel boost one might use a combination of two boost: along pY and along . 
We have checked numerically that both methods are fully equivalent. 
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emission can be expressed as a linear combination of the LO matrix elements 20 21 . The 
latter can be calculated using some effective incoming on-shell partons four-momenta, e.g. 
with the help the spin amplitudes. However, care must be taken while constructing these 
effective four-momenta for each LO matrix element individually in order not to spoil 
spin correlations. Coefficients of these matrix elements can be expressed as functions of 
variables related to the radiated partons (e.g. their momenta fractions and polar angles 
in an appropriate frame) and they are generally different for each LO matrix element. In 
a Monte Carlo approach, computations of the NLO cross section can be done with the use 
of the so-called branching algorithm, where in each branch a single LO matrix element 
is evaluated and a particular branch is picked up with a probability proportional to the 
coefficient of this matrix element. What is also interesting, flavours of these effective 
partons depend only on the respective LO process and can be different than flavours 
of partons initiating the NLO process. Recently, such a method has been adapted to 



implement the Drell-Yan processes in the POWHEG generator 22 . 

For the DY processes the differential cross section corresponding to the NLO real- 
parton radiation can be expressed as the following combination of the LO matrix elements: 

i 

where Ci are the coefficients depending on the radiated parton variables, e.g. Ci = 
Ci{x, cos 6), where x and 6 are the momentum fraction and the polar angle of this parton. 
The effective four-momenta Pq and Pq> of the incoming on-shell quark q and antiquark q' en- 
tering the LO matrix element are constructed in the NLO-process CM frame in such a way 
that the momentum of the "spectator" of the radiation is only rescaled without changing 
its direction, while the four-momentum of the "emitter" is calculated as a difference of the 
electroweak boson four-momentum and the "spectator" effective four-momentum, so its 
direction is different from that of the original "emitter" parton. The "emitter-spectator" 
assignment to the incoming partons is done on the Feynman-diagramatic basis, details 



can be found e.g. in 21 , 22j. 

In order to compute the appropriate LO matrix elements and to match the LO kine- 
matics with the NLO one, appropriate Lorentz transformations should be performed for 
the effective on-shell quarks four-momenta. Actually, they are analogous to the ones de- 



scribed above for the WINHAC-PYTHIA matching, see 22 . Based on this analogy we 
believe that our procedure for the kinematical matching between WINHAC and PYTHIA 
should be correct, at least up to the NLO. Of course, the parton shower provides only a 
leading-log (LL) type approximation of the NLO QCD corrections, but in PYTHIA the 
exact NLO matrix elements for real-parton emission can be taken into account through 
appropriate correcting weights. If this is done, the predictions of PYTHIA for the lep- 
tonic distributions in the DY processes should be exact at the NLO for the hard process, 
except for the normalisation. The latter will not be correct because PYTHIA does not 
include the NLO virtual corrections - this could be easily fixed by applying the NLO 
i^-factor. For the lepton charge asymmetries this i^-factor is even not needed because it 
cancels out between numerators and denominators. In the discussed WINHAC to PYTHIA 
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interface these NLO matrix element corrections are included by default. However, the 
PYTHIA predictions will not be exactly the same as from the fixed-order NLO calcula- 
tions, because PYTHIA generates also the higher-order LL-type QCD corrections through 
parton showers. They will lead to additional distortions of the leptonic distributions, in 
particular that of p^, but should not change them drastically. 



3 Numerical results and discussion 

An observable which is most sensitive to details of the kinematical matching between the 
ME-type MC generator and the PS generator in the charged-current DY process is the 
final-state charged lepton transverse momentum pip. This is because its distribution is a 
strongly varying function with a sharp Jacobian peak. Its shape is considerably affected 
by the non-zero VT-boson transverse momentum p^, see e.g. |5]. Moreover, since is a 
vector boson and its coupling to fermions are of the V — A type, the angular distributions 
of its decay leptons in the H^-rest frame are highly asymmetric: 

da ( \^ 

oc \\-Qw^o^Qiq] , (2) 



d cos 9 



Iq 

where 9iq is, in the limit of massless quarks, the angle between the outgoing charged lepton 
and incoming quark directions, and Qw = ±1 is the ly-boson electric charge in the units 
of the positron charge. Because of that, pip depends strongly not only on p^ but also on 
the individual momenta of the quark and antiquark. 

In our numerical comparisons we have used the following MC programs: WINHAC 1.35 



(12|, PYTHIA 6.401 [23], MCONLO 4.03 |24]| and MCFM 5.8 [25 . 

The lepton charge asymmetry observables are used to scrutinize the differences be- 
tween the and W~ mediated processes in these generators. For a given kinematical 
variable a, the charge asymmetry Asym (a) is defined as: 

da~^/da + da~ /da 

where -|- and — refer to the electric charge of the W boson (or the final-state charged 
lepton) and da^ /da is the differential cross section of an observable a. 

The asymmetry distributions have been obtained for the proton-proton collisions at 



s = 14TeV using the CTEQ 6.1 PDF set 26 and the particles properties from the 



PDG 2011 publication 27 , for the following two cases: (1) without any kinematical 



restrictions for the outgoing lepton and (2) with the kinematical cuts 

pip > 20 GeV, \7]i\ < 2.5, E^''' > 25 GeV. (4) 



In Fig. [l] we show the charge asymmetry distributions as a function of pip for electrons 
obtained from WINHAC interfaced with PYTHIA. The kinematical matching described the 
previous section was applied with the effective on-shell quarks four-momenta as provided 
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WINHAC vs. PYTHIA: no cuts WINHAC vs. I\/IC@NL0: no cuts 




p; [GeV] p; [GeV] 



Figure 1: The comparisons of the pip charge asyrametries from WINHAC and PYTHIA 6.4 (left plots), 
and from WINHAC and MCONLO (right plots); lower plots show the differences between the programs. 



by PYTHIA. These distributions are compared with the ones coming directly from PYTHIA 
6.4 (left plots) and the ones obtained from MCONLO (right plots). Lower plots show the 
differences of the distributions presented in the upper plots. A very good agreement 
between WINHAC and PYTHIA shows that all the technical aspects of the kinematical 
matching in their interface were done correctly. However, we see that the WINHAC and 
PYTHIA results differ considerably from the MCONLO results, in particular we observe 
the opposite behaviour of the asymmetry above the Jacobian peak (> 40 GeV). The 
region around the Jacobian peak is crucial for the W mass measurements at the LHC, 
see e.g. Refs. [l|[3]. The source of this discrepancy must be understood to hope for any 
improvement of the present precision of the the W-hoson mass, width and their charge 
asymmetries at the LHC. 

In Fig. [2] we present the comparisons of the p^j. charge asymmetry distributions for 
muons without kinematical cuts and with the cuts of eq. (|4]) between WINHAC (with the 
same kinematical matching as in Fig. [T| and MCONLO. Similar discrepancies as for the 
electrons are observed for the fully inclusive distributions. In the presence of cuts the 
asymmetry distribution changes its shape and the differences between the two programs 
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Figure 2: The p^j. charge asymmetries from WINHAC and MC@NLO without cuts (left plots) and with 
the typical ATLAS and CMS cuts (right plots); lower plots show the differences between the programs. 



are smaller but still unacceptable. 

We have also compared the charge asymmetry distributions as a function of: , yw 
and rji, where: 



w 
Pt 



yw 



In 



2 -p^ 
-In (tan(e,/2)) 



w 



(5) 

(6) 
(7) 



and found a good agreement between WINHAC and MCQNLO. In Fig. [3]we present charge 
asymmetries for rji. Except for the large values of rji, i.e. except for the region which is 
beyond the measurement domain of the ATLAS and CMS experiments, the agreement 
between the two programs is good. The large discrepancies between the two programs in 
the restricted phase-space, as specified by eq. (|4]), are thus important only for the charge 
asymmetries of the transverse lepton momentum distributions. 

One may argue that these discrepancies result from differences in the shape of the 
p^ distributions between WINHAC and MCONLO. Indeed, we have found that the 
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W* LI* v„: no cuts W* ii- v- with cuts 
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Figure 3: The rji charge asymmetries from WINHAC and MC@NLO without cuts (left plots) and with 
the typical ATLAS and CMS cuts (right plots); lower plots show the differences between the programs. 



distributions differ for WINHAC and MCONLO, but mainly at low (< 6GeV), while 
for higher values their ratio is flat. In the PYTHIA PS algorithm used by WINHAC the 

distribution at low values is affected mainly but the so-called intrinsic partonic fcr 
which is generated from a Gaussian distribution with an adjustable width. We have used 
this dependence and generated the samples of events with amplified differences of the 

distributions between the above two generators in the range which is well beyond 
the present measurement uncertainties. We have observed that the corresponding charge 
asymmetries of the pip distributions remained hardly changed. We have also compared 
these asymmetries for p^ > 6 GeV, where the ratio of the p^ distributions from WINHAC 
and MCONLO is flat, and found similar results. Finally, have checked that, in spite of 
differences in the absolute pj^ distributions, their charge asymmetries agree very well 
between the two programs. This proofs that the differences in p^ do not explain the 
large discrepancies in the pip asymmetries between WINHAC and MCONLO. Thus, the 
latter must be attributed to the differences in the effective polarisation of the WINHAC 
and MCONLO 14^-bosons. 

Can we find simple physical arguments to explain the shape of the pip charge asym- 
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Figure 4: The pip charge asymmetries from WINHAC with = and MC@NLO, without cuts (left 
plots) and with the typical ATLAS and CMS cuts (right plots) ; lower plots show the differences between 
the programs. 



metry distribution? Can we say which program is right and which is wrong? In order to 
try to answer the above two questions, we first produced the WINHAC distributions for 
the case of pj^ = 0, i.e. without the PYTHIA parton shower in WINHAC (using purely 
beam-collinear quarks from the standard PDFs). These distributions are shown in Fig. |4 
and compared with the previous MCONLO result. As one can expect, for the WINHAC 
= case the charge asymmetry distribution is flat below the Jacobian peak position 
40 GeV), and then rises very slowly with increasing pip. Its average value below the 
peak position reflects the difference between the total cross sections for positively and 
negatively charged DY processes which is driven by the effective excess of the w-quarks 
with respect to the rf-quarks producing the H^-bosons at the LHC. For the pip values 
above the Jacobian peak VT-bosons must be off-shell if p^ = 0. Since higher invariant 
mass prefers harder quarks and since u is on average harder than d, the relative number 
of produced rises with respect to W~ . In the following, this effect will be called the 
isospin effect. 

In the presence of the kinematical cuts we see a good agreement between the two 
programs below the Jacobian peak, and the discrepancy begins above this peak but is 
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smaller than without cuts and also than when PYTHIA is used. The sharp cut at = 
25 GeV for the WINHAC results comes from the cut on E^^'^'^ (in the case of = they 
are equivalent). It is rather striking that below the Jacobian peak the charge asymmetry 
of the pip distribution is generated mainly by the cut on the lepton pseudorapidity, rji, 
and is hardly sensitive to the W^-bosons transverse momentum spectrum. 

What is the reason for the observed shape of the distribution for = 0? In the 
W-rest frame the events with high pip correspond to r]i ~ 0, while the ones with low pip to 
large positive and large negative rji. If we take the +z axis along the quark momentum, 
then the W^-bosons will have preferably positive rapidity in such a frame, since the quarks 
are on average harder than the antiquarks. Thus, when we perform a boost to the LAB 
frame in the presence of symmetric cuts on rji, the events with negative rji will migrate 
in while the ones with positive rji will migrate out of the selected kinematical region (the 
VT-boson rapidity just adds to the lepton pseudorapdity). Since for W~ charged leptons 
are emitted preferably along its direction, then more events with low pip will move out 
than move in, while for it will be opposite. This is why we observe the decrease of 
the asymmetry with increasing pip up to the value close to the Jacobian peak position. 
Close to the peak position majority of leptons must have rji in the W-rest frame and 
the migration mechanism discussed above can be neglected. For such plj. values and above 
the discussed earlier mechanism related to the relative hardness of the distributions of the 
u and d quarks takes over and the asymmetry rises. 

Having understood the influence of the migration and the quark-isospin effects on the 
charge asymmetry distribution, let us try to answer our main question: do we understand 
the pip charge asymmetry when p^ > 0? 

As discussed before, the shape of the distribution below the Jacobian peak position 
is determined by the migration mechanism and is hardy dependent on the underlying 
distribution of p^ . Therefore, in this region, our previous analysis holds. We thus con- 
centrate on the region of the large pip (above the position of the Jacobian peak). The main 
difference here with respect to the p^ = case is that in addition to the isospin effect 
another effect come into play and become dominant: the effect of hard QCD radiation 
which influences the effective polarisation of the W^-bosons. In the discussion presented 
below the ly-polarisation is specified in the reference frame in which the spin quantisation 
axis is parallel to the direction of the W^-boson. 

It has been shown recently that for the processes of ly + jets production at the LHC 



left-handedly polarized Ws dominate over the right-handedly polarized l^s 28 . For 
the left-handed Ws the charged leptons are emitted preferably in the W~ direction and 
opposite to the W~^ direction (and vice versa for the right-handed Ws). Therefore, the 
non-zero p^ increases, on the average, the transverse momentum of the negatively charged 
lepton and decreases it for the positively charged one. This is what we observe in the left 
plot of Fig. E^where for MCONLO the asymmetry decreases for high pip. There is, of course. 



some contribution from longitudinally polarized Ws, but it never dominates 28 and, 
what is more important, charged lepton angular distributions are in this case identical 
for W~^ and W~. Moreover, the isospin effect, which could potentially counterbalance 
such a decrease, is sizeably smaller in magnitude due to a steeply falling Breit-Wigner 
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distribution. Therefore, the MCONLO results do have a rather convincing explanation of 
the charge asymmetry behaviour while the PYTHIA results do not. 
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Figure 5: The pip charge asymmetries from WINHAC and MCFM for W + I jet (left plots), and MCFM 
for W + 2 jets (right plots); lower plots show the differences between the programs. 



In order to check the validity of the above reasoning we have generated events using the 
Monte Carlo program MCFM |25] which calculates the fixed-order QCD corrections to the 
hard process convoluted with the coUinear PDFs. The charge asymmetry distribution for 
W + ljet and W + 2 jets are presented in Fig. |5] and compared with the ones from WINHAC 
with the standard PYTHIA parton-shower matching. The asymmetries predicted by 
MCFM are close to those from MCONLO, which supports the conclusion that the MCONLO 
predictions on the pip asymmetries are more likely to be correct than those of PYTHIA. 

But can we find the simple reason why the PYTHIA predictions are so grossly wrong? 
From our numerical tests and discussion presented above it becomes obvious that the 
problem must be related to the modelling of the effective polarisation of VT-bosons. In 
the LO approximation and for on-shell partons the VT-polarisation is uniquely driven by 
the asymmetry in the distributions of the momenta of the effective quark and antiquark 
entering the DY processes (see eq. ([T])), rather than by their sum which determines yw 
and . Inspecting the PYTHIA 6.4 manual we have found that the construction 
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Figure 6: The pip charge asymmetries from WINHAC with the transverse momenta of the effective 
quarks swapped and MCONLO, without cuts (left plots) and with the typical ATLAS and CMS cuts 
(right plots); lower plots show the differences between the programs. 



22 



of these effective on-shell partons momenta should agree at NLO with that of Ref. 
PYTHIA, of course, generates through the parton-shower more than a single NLO emis- 
sion, however they should not change considerably (or even revert) the NLO effective 
partons momenta as such additional emissions are mainly soft and collinear. Therefore, 
the PYTHIA predictions for the asymmetries should not differ much from the NLO 
ones. 

At this point we started searching not only for possible conceptual but also for the 
technical errors affecting the spatial orientation of the quark and antiquark momenta. 
We have made several technical checks of the PYTHIA generator along this line. One 
of the checks done was to swap the transverse momenta of the effective on-shell quark 
and antiquark. To our great surprise, once this was done on the event-by-event basis, we 
have obtained a very good agreement with the M CO NLO charge asymmetry distribution, 
both in the full phase-space and in the restricted kinematical region. The comparisons 
are shown in Fig. [g] for the pip and in Fig. [t] for the rji dependence of the lepton charge 
asymmetry. This agreement may be accidental but it may also suggest that the transverse 
momenta are, perhaps, not correctly assigned to the effective quark and antiquark in 
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Figure 7: The r/i charge asymmetries from WINHAC with the transverse momenta of the effective quarks 
swapped and MCONLO, without cuts (left plots) and with the typical ATLAS and CMS cuts (right plots); 
lower plots show the differences between the programs. 



PYTHIA. Whether or not such a hypothesis is true can however be verified only by the 
authors of the PYTHIA generator. 

On the conceptual side we have investigated the mechanism which drives the effec- 
tive LO polarisation of the W bosons in the DY process involving on-shell quarks. We 
have found that skipping the rotations of leptons momenta in the IV-boson rest frame, 
described in point 5 of the previous section, gives better agreement of the PYTHIA pip 
charge asymmetry with the MC@NLO one. The results are shown in Fig. |8| The agree- 
ment with MCONLO is slightly worse than in Fig. [6} however much better than in Fig. |2} 
Note that skipping these rotations is equivalent to retaining the PS-initial (parton-shower 
unaffected) effective on-shell quark helicities rather than those corresponding the PS-final 
ones (following the parton-shower). 

We have implemented the above two options in the new version of WINHAC 
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These versions cannot replace the future state-of-the-art NLO programs with the NLO 
PS and the full set of EW radiative corrections. However, as long as such programs are 
not available, they may be of use in the initial phase of the measurement of the lepton 
charge asymmetries at the LHC First of all, they can be of help in the unfolding of the 
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Figure 8: The p^p charge asymmetries from WINHAC without rotations of leptons momenta in the 
VF-rest frame and MCONLO, without cuts (left plots) and with the typical ATLAS and CMS cuts (right 
plots); lower plots show the differences between the programs. 



measured charge lepton asymmetries in the experimental procedures where the precision 
of the EW and the real photon radiative correction matters. More importantly, the above 
versions, providing the simplified LO picture of the effective polarisation of ly-bosons 
at the LHC, may help in designing the new polarisation-dependent observables for the 
studies of the electroweak symmetry breaking (EWSB) mechanism. 



4 Summary 

In this paper we have discussed the generic problem of kinematical matching of a parton 
shower generator with a matrix element generator for the Drell-Yan processes involving 
spin 1 intermediate particles. We have argued that the Les Houches Accord must be 
extended to take into account the spin correlations at all the stages of the event generation. 
We have described in detail our kinematical matching procedure which is used in the 
interface of our WINHAC generator with the PYTHIA 6.4 generator. We have demonstrated 
that the momentum vectors of the on-shell quark and antiquark, carrying in the LO 
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approximation the full information on the W^-boson polarisation, must be well defined. 
Any error in directions of these vectors has dramatic consequences for the dependence 
of the charge asymmetries at the LHC. In particular, using the transverse momenta of the 
effective quarks provided in PYTHIA 6.4 leads to completely different behaviour of the 
above asymmetry than predicted by the NLO (and beyond) calculations, e.g. MC@NLO 
and MCFM. We have found that simple swapping of the effective quark and antiquark 
transverse momenta in PYTHIA, or skipping the rotation of the outgoing lepton momenta, 
results in the pip charged asymmetries that match the NLO predictions of MC@NLO. We 
have implemented the corresponding matching schemes in the new version of the WINHAC 
generator. 

The issue of the proper matching between the matrix element calculations and the 
parton-shower generators respecting the spin correlations is important not only for the 
charged-current Drell-Yan processes but also for any process of production and decay of 
non-zero spin particles at the LHC. It needs to be readdressed in the more general context 
of matching the NLO matrix elements with the NLO parton shower such that a handle is 
given to the experimentalists to control the relative contributions of all the spin density 
matrix elements of the decaying particles, thus allowing for an experimental verification of 
the implemented Monte Carlo mechanism which drives the polarisation of non-zero spin 
particles at the LHC. 
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